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Abstract 

We present in this paper a simple method to obtain informations on 
galaxy models parameters using the measured value of the microlens- 
ing optical depth. Assuming a 100% MACHO's dark halo, we ask the 
predicted optical depth for a given model to be the same as the ob- 
served one, in a given direction. Writing the optical depth in terms 
of the given halo model parameters and inverting this relation with 
respect to one of them, it is possible to get information on it, fixing 
under reasonable hypothesis the other parameters. This is what we 
call the inverse problem in microlensing. We apply this technique to 
the class of power -law models with flat rotation curves, determining 
the range for the core radius R c compatible with the constraints on 
the halo flattening q and the measures of r towards LMC. Next, we 
apply the same method to a simple triaxial model, evaluating the axial 
ratios. 
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1 Introduction 



Gravitational microlensing is a powerful tool to investigate the dark halo, 
whose presence around galaxies is suggested by the flatness of rotation 
curves. The discovery of microlensing events towards LMC and SMC has 
witnessed the presence of MACHO's as baryonic costituents of the dark halo, 
but has also opened many new problems. What is the fraction / of dark 
matter composed by MACHO's? How are they distribuited? And (last but 
not least) what is the shape of the dark halo? 

In these years several authors have tried to answer these questions using 
the data coming from microlensing observations towards Magellanic Clouds 
to learn something about the properties of the dark halo. To start the inves- 
tigations, Paczynski (1986) has adopted the most simple model: the singular 
isothermal sphere. Few years later, Griest (1991) suggested to add a core 
radius to eliminate the central singularity of the mass density of Paczynski' 
s model, obtaining what is now called the standard model. Unfortunately, 
optical depth and event rate so predicted are higher than observed, therefore 
many authors have concluded that MACHO's are not the only components 
of the dark halo. A simple way to estimate / is to do the ratio between the 
observed optical depth and the predicted value for a model totally made of 
MACHO's. However, this method makes / model dependent, i.e. / depends 
on the model adopted to describe the shape of the halo and the distribution 
of MACHO's inside it. In principle, there is no reason a priori to say that 
MACHO's are not the only components of the dark halo. In this paper we 
consider / = 1 and show how it's possible to use the measured value of 
the optical depth to extract informations on galaxy parameters. Instead of 
going from galactic dynamics to microlensing, we try the other way around, 
that is going from microlensing to galactic dynamics. We will express the 
optical depth r as function of some galaxy model parameters (x,y, ...) and 
then invert the relation so found. If / = 1, then Ttheor, he. the value pre- 
dicted for a 100% MACHO's halo, has to be equal to the measured value 
Tabs', so, we may use the function r = r(x,y, ...) to extract informations on 
(x,y, ...). This is what we call the inverse problem in microlensing^ 

In this paper, we apply this method of investigation to two classes of 
galaxy models and try to learn something about their characteristics. We 
do not consider spherical models since there are few evidences which lead to 

X A similar technique has been yet used in other works ( see e.g. De Paolis, Ingrosso & 
Jetzer, 1996), but it has never been applied to the models we consider here. 
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reject this extremely simple shape for the dark halo. In fact, it is possible 
to estimate the halo flattening through dynamical modelling of collissionless 
kinematic tracers, such as extreme Poulation II stars. These probe the 
potential at large distances and high above the disk and permit to estimate 
the ratio c/a between the principal semiaxes, which turn out to be (c/a)t rc = 
0.6 - 0.85 (| ommer-Larsen Zen, 1990|) . Using the Jeans equations, van 
der Marel (1991) has shown that halo flattenings of 0.4 < (c/a)haio 

< 1 

were consistent with the data depending on the tilt of the velocity ellipsoid 
away from the disk plane. These evidences suggest to give off the spherical 
simmetry. A step in this direction is the adoption of axisimmetric model for 
the dark halo, this is the reason why many authors have used the so called 
power -law models ( [Evans, 1993| ), ( Evans, 1994]) q Optical depth and event 
rate towards LMC have been calculated ( [Evans fc Jijina, 1994| ), (Alcock 



et al., 1995) and they found a method to experimentally estimate the halo 



flattening from the ratio between event rate observed towards LMC and 
SMC. Unfortunately, predicted and observed optical depth are once again 
not consistent. In this paper, we apply the inverse problem method to 
the class of power - law models with flat rotation curves and try to learn 
something about their properties to see whether it is possible to explain the 
disagreement yet found between theory and observations of microlensing 
events towards LMC. 

It is worthy to note, however, that there are no direct evidences of de- 
viations of the Milky Way halo from axisymmetry in the disk plane. This 
may be the consequence of a partial compensation between disk and halo 
potentials, in case they are of comparable dynamical importance. It happens 
infact that, even for significantly triaxial halos, the combined disk -halo po- 
tential is much more symmetric than the single components potentials ( [Rix 



1996] ). Then, it is not clear whether one must give off the axisymmetric 



hypothesis. An help in this choice may come from the theoretical predic- 
tions for the halos arising in hierarchical structure formation scenarios. For 
non dissipational dark matter, this question has been addressed by cosmo- 
logical N-body simulations ( pubinski fc Carlberg, 1991 ), (Katz &: Gunn 



1992| ), ( [Navarro et al., 1996| ). Despite the ongoing discussion about numer- 



ical simulations and initial conditions, all dissipationless simulations agree 
on predicting a strongly triaxial shapes, closer to prolate than oblate. The 
resulting axis ratios vary slowly with radius and the distribution of halo 



H owever, flattened models, other than power -law, has been also considered (Jetzer 
199l|). 
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flattenings (c/a) peaks near 0.7. However, dissipational dark matter (such 
as MACHO 's and cold molecolar clouds) will affect the halo shape in dif- 
ferent ways, compressing the halo in the disk plane direction. Simulations 
including gas dynamics ( Katz fc Gunn, 1992|) or a slow growth of a disk 
component ( Dubinski, 1994] ) have found that halos become slightly flatter 
and much closer to oblate and that these changes extend to large radii. Our 
position in the Milky Way complicates measurements of the halo shape so 
that many models are consistent with the data coming from different kind of 
measurements (for a brief review see ( Rix, 1996 )). Better results have been 
obtained for external disk galaxies. In particular, flattening estimates are 
possible for several polar ring galaxies ( see e.g. ( Sackett fc Sparkle, 1990 )), 
whose halos are definetely flattened perpendicular to the stellar disk accord- 
ing to numerical simulations. If our dark halo is typical, then cosmological 
simulations and observations around other spiral galaxies suggest a triaxial 
halo, flattened perpendicular to the disk plane. A result in this sense has 
been obtained by Binney, May & Ostriker (1987). They constructed approx- 
imately self - consistent disk halo models with Staeckel potentials and find 
that (c/a) halo ~ 0.5 fits the observations best. 

Following the suggestions coming from numerical simulations and the 
cited result of Binney, May &; Ostriker, in this paper we also adopt a par- 
ticular separable triaxial model to describe the dark halo and evaluate the 
predicted optical depth towards LMC. Doing this requires the knowledge of 
the axial ratios and the extension D = a\ of the model, which are unknown. 
However, we will show how our method may estimate the axial ratios under 
the hypotheses / = 1 and D = 100 kpc. 

In Section 2, we briefly review some basics notions of microlensing, intro- 
ducing the optical depth r and the event rate V. Section 3 is devoted to the 
description of power - law models, while in Section 4 we calculate the optical 
depth towards LMC and apply the method of the inverse problem. Section 
5 is dedicated to the discussion of the results obtained for the power -law 
models with flat rotation curves. The perfect ellipsoid, a particularly simple 
separable triaxial model, is described in Section 6, while in the following 
Section 7 we do the same analysis as in Section 4 for this model. The rela- 
tive results are summarised in Section 8, while the next section is devoted 
to conclusions. 
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2 Basics of microlensing 



When a lens moves across the line of sight, the light coming from the source 
is amplified by a factor ( [Morlerach fc Roulet, 1997| ), flJetzer, 199S| ): 

A(u) 



uVu 2 + 4 

where u = b/R e is the impact parameter in units of Einstein radius : 



R P 



UGMD d (D s - D d ) 



Here, D s is the distance of the source, is the distance along the line 
of sight and M is the mass of the lens. If the MACHO is moving with a 
transverse velocity v±, the impact parameter is time dependent and is given 
by: 



u(t) = ^u 2 nm + u 2 (t-t ) 2 (3) 

being u m i n = u{A max ) and ui = v±/R e . Experimentally microlensing events 
are characterized by the maximum amplification A max and the event dura- 
tion T = 1/uj. 

Two are the most important parameters in microlensing analyses : the 
optical depth and the event rate. According to Vietri &: Ostriker (1983), 
the optical depth r is the probability of observing a microlensing event. 
In the little number limit, this is just the number of MACHO's inside the 
microlensing tube, a cylinder whose axis is the line of sight with radius 



UfhR e , being uth = u{Ath) with Ath the threshold amplification (Paczynski 



1986| ). When A t ^ = 1.34, = 1. The optical depth is then (Paczynski 



T9861) , flJetzer, 19981) : 



r = *J jjPRlds (4) 

where we have posed s = Dd/D s . 

The event rate V is simply the number of events divided by the duration 
time of observation. Analitically, it is defined as ( De Rujula et al, 1991| ) : 



e(T)§*T (5) 
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where s(T) is the efficiency function and dT/dT is the differential rate; in 
the hypothesis of same mass, it is given by : 



dT r 2n r Ds / R \ 4 r 



being /(x, v j the distribution function^ of the model adopted for describing 
the distribution of MACHO's and vn is the MACHO velocity parallel to the 
line of sight. 

Optical depth, event duration and event rate are related by the useful 
relation (Dc Rujula et al, 1991), ( Jetzer, 1998| ) : 

<r>=-^ (7) 

where we have defined : 

<T> =fJ T § dT - < 8 > 

Eq. ([?]) is similar to the one used to estimate optical depth from the observed 
events. The latter is : 

- _ V Ti (Q\ 
Tobs ~ 2iV*U* <T t ) W 

where the sum is on the events observed, being iV* the number of resolved 
monitored stars, t Q b s the duration of observations and e(Tj) the efficiency for 
the i- th event. 



3 Power - law models 

The event rate is not easy to calculate since the DF /(x, v) cannot be pre- 
scribed arbitrarly. Under reasonable hypotheses, the halo is a collisionless 
system and its DF must satisfy collisionless Boltzmann equation ([Binney 



fc Tremaine, 198^ ). By Jeans' theorems, the DF depends on configuration 



space coordinates only through the three (or just one or two of them) iso- 
lating integrals of motion in the gravitational potential generated by the 
galaxy itself. So the construction of self consistent galactic models is not so 
easy, especially for non spherical systems. 

3 Remember that the distribution function (hereafter DF) of a collisionless system is the 
numerical density on the configuration space, i.e. the number of particles with position in 



the interval (x, x + d x) and velocity in (v, v + d v) (Binney & Tremaine, 1987) 
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We start to dedicate our attention to a class of axisymmetric galactic 



halo models: the power -law models, also known as Evans models ( Evans 



1993| ) ( [Evans, 1994 ). In this case the two integrals of motion are the energy 



E and the z - component L z of angular momentum. The DF is : 

f(E, L z ) = + BB A/P~l/2 + CE 2/P-l/2 (for > Q) . 

f(E,L z ) = AL 2 z exp4E/vo+Bexp4E/vo + Cexp2E/v% (for (3 = 0) ; 
AL 2 B C 

f^ L ^ = (-3)3/2-4/, + (_ g) l/2-4/, + { _ E) l/2-2 /0 (for/?<0)(10) 

where A, B, C are constants given in ( [Evans, l~994l) . 

Integrating the DF on the velocity space, one obtains the mass density 
of the models, which is all we need to evaluate the optical depth : 

m v *gig R 2 c (i + V) + R 2 (i - Pq 2 ) +z*[2- g -2(i + /?)] 



where (i?, z) are the usual cylindrical coordinates. The rotation curve is 
given by : 

""t ni 

The parameters of the power - law models are the following : 



„,2 n/3 r>2 

vt(R) = ? C , +2 (12) 



1. The core radius R c , wich measures the scale at which the density 
begins to soften. Its value is not well known; Bachall, Schmidt & 
Soneira (1983) give an estimate of R c as 2 kpc from star count data, 
while Caldwell & Ostriker (1981) state R c = 10 kpc. However, values 
as large as 20 kpc are also possible if we accept the hypothesis of 
maximal disk, i.e. rising the total mass of the disk and its contribution 
to the gravitational potential of the whole Galaxy. 

2. The parameter f3, which determines the asymptotic behaviour of the 
rotation curve; in fact, it is defined as : 

dlog^QB) 

(5 = -lmiR^oo — (13) 

The value of f3 is not known since it is not possible to measure the star 
rotation curves of the Galaxy beyond 20 kpc. However, the data seem 
to suggest that Milky Way is similar to other spiral galaxies, whose 
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rotation curve is flat {(3 = 0) or only slowly rising(/3 < 0). So, in the 
following we confine our attention to the case (3 = and try to learn 
something on the other parameters from the analyses of the present 
measure of optical depth ( see later, Section 4). 

3. q describes the flattening of the dark halo : actually it is the axial 
ratio of the equipotential spheroids; q = 1 is the value for a spherical 
halo, while q = 0.7 is for an E6 halo flEvans, T994p . The isophotal 
ellipticity is a function of q and of the other parameters of the model. 
Obviously, no visible constraints are known for the dark halo; so we 
have to examine the whole physically possible range for q. This is fixed 
by the requirement that the mass density and the DF are always not 
negative. If we choose (3 = 0, these conditions define the range flEvans, 
1993|) : 

0.71 < q < 1.08. 



4. The distance Rq of the Sun from the Galactic centre. Recently, this 
quantity has been reviewed by Reid (1989), who found that most recent 
determinations lie between 7 and 9 kpc; we will use the value 8 kpc, 
but we also consider the other values. 

5. The normalization velocity v$ which determines the typical velocity of 
MACHO's in the halo. This value is fixed by the local circular velocity, 
v c (Rq), through the formula (|i~2|). In our case ((3 = 0) we have : 



Rn 



< = -^r^vi(Ro) (14) 



We have to choose a value for v c (Rq); following Merrifield (1992), we 



adopt v c (Rq) = 200±10 km/s, even if the IAU value is 220 km/s ([Kerr 
fc Lynden- Bell, 1986Q . 



4 The optical depth for power - law models with 
flat rotation curves 

To estimate the optical depth as in eq. @), we need just the mass density 
of the models, which is given in (111]). Changing coordinates from cylindri- 



cal to heliocentric and executing the integral, one obtains (Evans fc Jijina 



1994) , ( |Alcock et al, 1995Q : 

r = v%R%< f 1 s(l-s)(A T s 2 + B T s + C T ) 2 dc 
<?q*Dl Jo ( DtS 2 + EtS + Ft) ^ 



where the constants A T , B T , C T , D T ,E T , F T are given in eq. A7 of (Alcock et 



al., 1995). For the class of power -law models with /3 = 0, they are : 

A T = cos 2 b + (2 - q~ 2 ) sin 2 b; 

B T = —2D~ 1 cos b cos I; 

C T = [R 2 + R 2 c (l + 2q- 2 )]/D 2 s ; 



D T = cos 2 b + q 2 sin 2 b; 



E T = —2RqD s 1 cos b cos/; 



F T = (i2§ + i£)/D|. 

The optical depth is then easy to evaluate in this case (in which (3 = 
and u t h = 1) and the result is : 

r = 4k (16) 
c z q z 

being I the integral in eq. (|l5|). 

In fig. la we report the optical depth towards LMCQ in function of the 
core radius R c and the flattening parameter q, having fixed : j3 = 0, Rq = 8 
kpc, v c (R ) = 200 km/s. 

Following the idea described in the introduction, we use this graph to 
extract information about the parameters of this class of power -law models, 
trying to learn something about the shape of the halo or the asymptotic 
behaviour of the rotation curve. The measured value of the optical depth 
towards LMC is : T obs = 2.1+^| x 10" 7 QBennett, 1998| ). We consider then 
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Figure 1: a.) Optical depth r vs. q and R c for power -law models with 
j3 = 0, R = 8kpc, v c {R ) = 200km/s; b.) section of a) for r = 2.1 x 1(T 7 . 

a section of fig. la at the level r = 2.1 x 10~ 7 and obtain the graph in fig. 
lb. 

For every value of q in the physically possible range, there is only one 
value of the core radius R c for which the model ( with the specified value 
of q and the other parameters fixed as before ) reproduces an optical depth 
equal to the one observed. In particular, the range for the core radius turns 
out to be : 

14.03 kpc <R C < 15.91 kpc (17) 

This may be compared to the range delimited by observations, which is 
2 -J- 10 kpc; excluding the hypothesis of maximal disk, the two ranges do not 
overlap. 

4 D B = 50 kpc; (b,l) = (-32° .9, 280°.5). 
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Figure 2: Optical depth r (towards LMC) vs. halo flattening q, having fixed 
P = 0, R c = 5.6 kpc; R = 8 kpc, u c (i2o) = 190, 200, 210 km/s. 

However, one must also consider different choices for v c (Rq) and Ro 
other than 200 km/s and 8.0 kpc. In fact, r is strongly dependent on this 
two parameters. In fig. 2, we report the graph of r vs. q for v c (Ro) = 
190, 200, 210 km/s, having fixed /3 = 0, R c = 5.6 kpc ( [letzer, 1998|) , R = 8.0 
kpc. 

The relative variation of r (for q = 0.8) is : 

. / x t(v c (Ro) = 210) - tCvJRq) = 190) 

Ar v (q = 0.8 = 1 Cl ° J , ' 1 ^ 20%. 

r(v c (i<o) = 210) 

Fig. 3 is, instead, a graph of r vs. q for R = 7, 8, 9 kpc, being v c (Rq) = 
200 km/ s and (3 and R c as before. The relative variation is : 

ATrh = . 8) = ^ = T)-,W, = 9) : 
r(i? = 7) 
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Figure 3: Optical depth r (towards LMC) vs. halo flattening q, having fixed 
P = 0, R c = 5.6 kpc; v c (R ) = 200 km/s, R = 7, 8, 9 kpc. 

Being the relative variations so strong, we do the same analysis as the 
previous one for every values of Rq and v c (Rq), evaluating in each case the 
range for the core radius R c . The results are summarized in table 1. 

It is worthy to note that the range for R c never overlaps the observational 
range and, in one case, exceeds also the 20 kpc permitted by the maximal 
disk hypothesis. 

Finally, we consider the observational range of T bsi which is, as we have 
already pointed out, 1.3 x 10 7 < T j, s < 3.4 x 10~ 7 ( |Bennett, 1998j ), and we 



do the same analysis for the minimum and the maximum value compatible 
with errors on T Q b s , limiting ourselves to the models with f3 = 0, Rq = 8 kpc, 
v c(Rq) = 200 km/s. For r equal to 1.3 x 10 we obtain : 

9.05 kpc < R c < 10.51 kpc , (18) 
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Table 1: Range for the core radius R c for models with (3 = and different 
values of Rq and v c (Rq). 



Ro ( kpc ) 


v c (R ) ( km/s ) 


R c ( kpc ) 


7 


190 


12.48 -7- 13.86 


7 


200 


11.37 -r- 12.83 


7 


210 


10.51 -T- 11.62 


8 


190 


15.57 -j- 17.55 


8 


200 


14.03 -7- 15.91 


8 


210 


13.00 -=- 14.63 


9 


190 


18.97 -7- 21.91 


9 


200 


17.24 -7- 19.91 


9 


210 


15.83 -T- 18.15 



which is compatible with the range outlined by observations, while none of 
our models is compatible with a value of r equal to 3.4 x 10~ 7 . 

5 Discussion of results for the power - law models 
with = 

We fix now our attention on the results obtained for different power - law 
models with (3 = and Rq and v c (Ro) compatible with observations, sum- 
marized in table 1 . The ranges for the core radius R c compatible with the 
hypothesis of a dark halo entirely composed of MACHO's are given in the 
third column of the table and we see that they never overlap the observa- 
tional range of 2 4- 10 kpc. We can draw some possible conclusions : 

1. If we want to keep the hypothesis of axisymmetry, we could take into 
consideration values for (3 different from zero, i.e. rotation curves 
which are no more asymptotically flat. In this case, one has to fix a 
criterium to consider either power -law models with (3 > (declining 
rotation curve) or (3 < (rising rotation curve) and then repeat the 
analysis present in this paper and then study the results. 

2. The halo is axisymmetric with flat rotation curves, but the disk is 
maximal, so that the core radius may be higher than 10 kpc. However, 
in this case also the contribution to the optical depth of the disk lenses 
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has to be taken into consideration and the extension of the dark halo[] 
has to be evaluated. 

3. The dark halo is not entirely composed of MACHO's. The same anal- 
ysis may be repeated making sections of fig. 1 at different values of r, 
choosing the fraction / of the halo in form of MACHO's and imposing 
that the predicted optical depth r for a halo with the specified / is 
equal to the measured value. 

At the moment it is not easy to choose between the possibilities outlined 
above. An important point is that, to investigate the shape of the dark halo 
it is necessary to make observations towards directions other than LMC; 
in this way more informations on the MACHO's distribution in the Galaxy 
would be at our disposal. 



6 Separable triaxial models and perfect ellipsoid 

Now, we turn our attention to another class of halo models, which is more 
general than the axially simmetric we have previously considered, the triaxial 
halos. Among the triaxial models, the simplest ones belong to the class 
of the separable triaxial models, a class of self- gravitating systems whose 
gravitational potential is of Staeckel form QStaeckel, 1890 ), ( ptaeckel, 1893| ). 



This kind of potentials are easier to study if we adopt ellipsoidal coordinates, 
which are defined as follows ( |Morse Feshbach, 1953] ) : 

x 2 _ + a )(/ i + at)(y + a) 



y 2 



2 

Z 



(a - 0){a - 7) 

(A + (3)(fi + P)(is + P) 

(P-a)(P -7) 
(A + 7)^ + 7)^ + 7) 



(7-a)( 7 -/3) 

where a, (3, 7 are three arbitrary constants. If two of them are equals, ellip- 
soidal coordinates degenerate in oblate (or prolate) spheroidal coordinates 



e Zeeuw, 1985a|) , while they reduce to the usual spherical coordinates 



when a = f3 = 7. 



5 Remember that in this paper we have taken an halo extended till LMC, which is a 
reasonable hypothesis. 
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In the case of a Staeckel potential there are three indipendent integrals of 



motion ( De Zeeuw, 1985a ) and these models have been extensively studied 
in literature. In a classical paper, Kuzmin has studied axysimmetric oblate 
systems with Staeckel potential and has demonstrated a very useful theorem 
flKuzmin, 1956| ), later generalized to generic separable triaxial models by 
De Zeeuw (1985b). According to this theorem, assigning the mass density 
profile along the z - axis is sufficient to determine the whole density profile 
of a family of models, all having, along the z-axis, the given mass density 
profile. Therefore, it is a relatively simple work to construct various models 



of separable triaxial systems with a Staeckel potential (De Zeeuw et al. 
"19861 ), ( Pe Zeeuw fc Pfenninger,"T988|) . 



We limit ourselves to a class of separable models, called perfect ellipsoids 
)e Zeeuw, 1985a|) , ( pe Zeeuw et al, 19861 ), (|De Zeeuw fc Pfenninger, 1988|) . 



We assume the following form for the density profile along the z - axis 



p c n 



p(0,0,z)=^)= ( /; c2)n (19) 

being c an arbitrary constant, po the central density, we consider only the 
case n > 0. Following ( De Zeeuw, 1985b| ), one may demonstrate that the 



mass density (in ellipsoidal coordinates) is : 

p(X, fi, v) = g\^{\) + g^ipi/J-) + 9u 2 ^iy) + 

{*(A) - *(m)} , {*(jj) -*!/')} 

" " T 

(20) 



+ 2g\9fi 7T x \-2gfj,g v r h 

(A - n) p \Ji-v) 



{v&H - tt(A)} 

+ 2 ^ gA („ - A) 
where (£ = \,/J,,v) is defined as : 



l 

i\ jt-i 







while we have posed : 



(A-mKA-i/) 
(/i + q)(// + /?) 

(/X -!/)(// -A) 



^ = f + t^ (22) 
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9v 



{y + a)(y + ff) 
(z/- A)(i/-/x)" 



Without loss of generality, we may choose c 2 
becomes : 

Poc n 



£n/2 



while its primitive is : 



*(0 



„n-2 



f(n-2)/2 



(23) 

-7 in (H), so eq. © 
(24) 



for ji / 2, 



\t(0=A)C 2 ln- 



for n = 2. 



The case of perfect ellipsoid corresponds to choosing n = 4 in 

(U); after 

some calculations, one gets the following expression for the mass density in 
ellipsoidal coordinates : 



(25) 



In (^) a, /3, 7 are arbitrary constants. Choosing 



= -a\ 7 
and adopting cartesian coordinates, eq. (|25| ) becomes : 

Po _ Po(a! + -Ro) 2 



«3 



(1 + m 



a|(l + fh> 



-n 2 \ 2 



(26) 



Here, p & and Rq are respectively the local mass density and the distance 
of the Sun to the galactic centre and : 



m 2 



x 2 y 2 z 2 

~ 2 3 72 " 



a 



Co 



(27) 



In ( |27|) ai,a2 and 03 are the semiaxes of the concentric equipotential 
ellipsoidal surfaces. They must satisfy the relation a\ > a 2 > 03, but their 
exact values are not known. It is interesting to note that this is the only 
separable triaxial model whose equipotential surfaces are exactly similar to 
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concentric ellipsoids ( De Zeeuw fc Lynden- Bell, 1985| ). The mass inside the 
ellipsoidal radius is : 



/ 777 \ 

M e u = 2-7T/9oaia2a3 arctanm I . (28) 

V 1 + 771/ 

The total mass is finite and is given by : 

M t ot = it 2 poa^az = it 2 p Q aia 2 a 3 (l + R\ja 2 2 ) 2 . (29) 
The DF of this model is not known. Some authors ( Dejonghe Sz Laurent] 



1991| ) have recently tried to reconstruct the DF of separable triaxial models, 



looking for a formula similar to the Eddington one for spherical models 
( Binney & Tremaine, 1987). The result is very difficult to obtain and it 



is not unique, this is the reason why we cannot calculate the event rate, 
therefore we limit ourselves to evaluate only the optical depth. We do it in 
the following section. 

7 Optical depth towards LMC and estimate of ax- 
ial ratios 

Evaluating optical depth through the formula (||) requires the mass density 
of the model, which is given, for the perfect ellipsoid, by equation (|26|). 
However, we have first to change the coordinates, from galacto centric to 
heliocentric and then use the formula (|j). The result is : 

4nGD 2 4 4/2, n2\2, fo n \ 

r - — 12 — aia 3 {a 2 + R )I (30) 



c- 



where I is the integral : 



1 = t FA ^ ^2 ds ( 31 ) 

Jo (As 2 + Bs + C) 2 y ' 



This may be calculated analitically. We get : 

R 4- 941 

(32) 



2 2B + AC 
I = -r + 



A A 3 / 2 



b B + 2A 

arctan —= — arctan 



A VA 



with 



A = D 2 (a2a 2 sin 2 b cos 2 I + of a 2 sin 2 6 sin 2 I + a 2 a 2 cos 2 b); 
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Table 2: Values of 02/01 and 03/01 to have a perfect ellipsoidal dark halo 
which predicts an optical depth towards LMC equal to what is measured. 



a 2 /a\ 


a-i/ai 


0.2 


0.064 


0.3 


0.063 


0.4 


0.062 


0.5 


0.061 


0.6 


0.061 


0.7 


0.061 


0.8 


0.061 


0.9 


0.061 


1.0 


0.061 



B = 2D s RQa\a\ sin 6 sin I; 

C2 2 r>2 i 2 2 2 
= a l a 3 K + a 1 a 2 a 3 ; 

A = AAC - B 2 > 0. 

Now, to calculate the optical depth towards LMC, we have only to insert 
in (g(]) and © the values (b,l) = (-32°.9, 280°.5) and D s = 50kpc; we 
get the expression of r in terms of the semiaxes of the perfect ellipsoid. We 
apply the method of the inverse problem to estimate the axial ratios of a 
perfect ellipsoidal dark halo for which the predicted optical depth is equal to 
the observed one. However, we have to fix the value of the halo extension, 
i.e. we have to fix D = a±. This latter quantity is not well determined and 
different authors give different estimates of D, in the range (50 -j- 200) kpc 
(for a brief review, see ( Zaritsky, 1998] )). A possible choice is D = 100 kpc. 



In fig. 4a we plot the optical depth towards LMC as function of axial ratios 
(Pz/ai, a 3 /a\, while fig. 4b is a section of 4a at the level r = 2.1 x 10~ 7 . 

Prom fig. 4b one can see that for each value of a 2 / 'a\ there is only one 
value of 03/01 for which the predicted optical depth equals the observed 
one. Table 2 shows how 03/01 changes varying 02/01. From the condition 
on ai, «2j 03 it has to be 03/01 < 02/01! from fig. 4b this happens when 
02/01 > 0.1, and the value of 03/01 assumes an approximately constant 
value. A good value is then 03/01 ~ constant = 0.06. 
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Figure 4: a) Optical depth r towards LMC vs. 0.2 1 a\ and a^/a\ for 
the perfect ellipsoid, having fixed the halo extension as a\ = D = 100 kpc. 
On the z-axis there is the quantity r r id = t/Gpq. b) Section of a) for 
r = 2.1 x 10~ 7 . 
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Now we should determine the other axial ratio, 02/01- We could use the 



formula (29) which gives the total mass of the halo, rewritten as : 



M „3/„ /„ x Cgo/gl + a 2/ Q l) /oq\ 

M tot = it p ai(a 3 /aij ■> — - — -7 (33 j 

(a 2 /ai) 4 

and determine the ratio 02/01 fixing the total mass of the halo. 

This latter is known with great uncertainties and different estimates are 
present, which adopt different values for the halo extension and a great 
variety of method and models (see e.g. ( [Lin et al, 1995[ ), ( [Zaritsky, 1998| )). 
We could use one of the estimates to obtain the corresponding value of 02/01, 
but we prefer just to point out the possibility of estimating the total mass 
of the dark halo, once the ratio 02/01 is known through other ways. Fig. 5 
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is a plot of M to t = Mtot((i2 / ai) , and only for 02/01 > 0.1 have to be taken 
into consideration. 



8 Discussion of results for the perfect ellipsoid 

Assuming a 100% MACHO's halo, we have analysed what are the features of 
a perfect ellipsoidal dark halo to have the theoretical optical depth towards 
LMC compatible with the observational values. We have obtained an esti- 
mate for the ratio 03 ja\ , and expressing the total mass of the dark halo as 
a function of the axial ratios, we have the possibility, either to estimate he 
mass of the dark halo, fixing a value for 02/01, or the axial ratio 02/01 once 
the total mass is known with a good degree of accuracy. The halo extension 
is taken as D = 100 kpc. 

The value we obtained for the axial ratio 03/01 ~ 0.06 shows an ex- 
tremely flattened dark halo in the direction perpendicular to the disk plane. 
This result is consistent with the predictions of numerical simulations with 
dissipational dark matter (Dubinski, 1994), ( [Rix, 1996| ). Ostriker & Peebles 



(1973) have studied the stability of flattened systems, determining a critical 
value for the ratio between kinetic and potential energy, above which the 
system is unstable against bar formation. They studied different models, 
but all of them are axisymmetric. This is not the case we are considering, so 
their results should be generalized in order to be applied. Few evidences are 



in favour of the presence of a bar in the bulge of our Galaxy ( Blitz fc Spergel 



1991|), (pwek et al., 199 ), some come from microlensing observations ( [Zhao 



et al, 1995[) , flZhao et al, 1996[) , but it is relative to the bulge component, 
and it is not easy to evaluate whether there is a significant contribution 
due to the halo, which is usually considered negligible in that region of the 
Galaxy. 



9 Conclusions 

In this paper we have seen how it is possible to infer the properties of the dark 
halo of our Galaxy from microlensing observations: what we have called the 
inverse problem in microlensing. This method of investigation of the dark 
halo has been applied to two different halo models, the power -law models 
with flat rotation curves and the perfect ellipsoid model. 

First, we have shown how the observed optical depth may restrict the 
range for the core radius for power -law models with flat rotation curves, 
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having supposed a 100% MACHO 's halo. The values we found are always 
higher than the observational ones, which are in the range (2 -j- 10) kpc, 
as delineated by the work of Bachall, Schmidt & Soneira and Caldwell & 
Ostriker. These results may be interpreted in different ways, as we pointed 
out in Section 5, and it is worthy to note that, when added to other ones 
coming from astronomical observations, our results give some indications on 
the validity of the description of the dark halo through power - law models. 

Next, we have considered a triaxial shape for the dark halo, adopting the 
perfect ellipsoid model, a simple separable triaxial system. Describing the 
halo as a triaxial structure is a difficult task since there are not enough 
observational constraints in order to determine all the properties of the 
model. However, fixing the halo extension, we have been able to determine 
one of the axial ratios (0.3/0,1) and find a relation between the total mass of 
the halo and the other axial ratio (02/01). 

We have applied the method of the inverse problem to constraint the 
dark halo on the basis of microlensing observations only in one direction, 
the LMC, and in fact we could not determine uniquely the properties of the 
models we considered. In order to investigate the shape of the halo and 
put more precise constraints on the relative parameters, it will be necessary 
to have microlensing data relative to observations towards other directions, 
like globular clusters, for example; they may be used as sources (see e.g. 
flGyuck fc Holder, 1997) )) or as sites of lenses when observing towards SMC 
(see e.g. flJetzer et al, 1998| )). The same method can be used, of course, to 
investigate other halo shapes. 

It is a great pleasure to thank G. Busarello and Ph. Jetzer for the 
discussion we had on the manuscript. 
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